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We have used optical tweezers to identify the elementary events underlying force generation in neuronal 
lamellipodia. When an optically trapped bead seals on the lamellipodium membrane, Brownian fluctuations 
decrease revealing the underlying elementary events. The distribution of bead velocities has long tails with 
frequent large positive and negative values associated to forward and backward jumps occurring in 0.1- 
0.2 ms with varying amplitudes up to 20 nm. Jump frequency and amplitude are reduced when actin 
turnover is slowed down by the addition of 25 nM Jasplakinolide. When myosin II is inhibited by the 
addition of 20 \iM Blebbistatin, jump frequency is reduced but to a lesser extent than by Jasplainolide. These 
jumps constitute the elementary events underlying force generation. 

Force generation is a fundamental process at the basis of cell motility 1 allowing neurons to explore the 
environment. Neuronal growth cones are the major motile structures located at the neurite tips 2 and are 
composed of lamellipodia and filopodia 3 . Lamellipodia are extended structures, from which filopodia 
emerge with a finger-like shape 4 . Their motion is essential during morphogenesis and for neuronal differentiation 
when their exploratory motion allows neurons to find the appropriate synaptic connections. Force generation is 
thought to be originating from the progressive addition of actin molecules to the existing network of actin 
filaments 5 and to be determined by the balance between actin polymerization and depolymerisation, modulated 
by controlling proteins 6 and by chemical and mechanical receptors coupled to the cytoskeleton 6-8 . However, very 
little is known about the elementary events underlying force generation. 

Actin polymerization has been primarily investigated in vitro by analysing the rate of elongation of isolated 
actin filaments. These investigations were performed with a low time resolution, often of the order of some tens of 
seconds and with a sensitivity of some hundreds of nm, providing values for actin polymerization rate ranging 
between 11.6 and 38 (1/uM s) 9 " 12 . 

Previous investigations in vivo using Atomic Force Microscopy 13 and opposing liquid flow 14 were limited to a 
temporal resolution in the 100 ms range and sensitivity of 50-100 pN. These experimental limitations can be 
overcome by using optical tweezers 15,16 , providing a ms resolution and pN sensitivity. In order to detect small 
displacements in the order of 2-5 nm it is necessary to reduce all perturbations by minimizing mechanical 
vibrations and performing the experiments under remote conditions (see Methods). By using these procedures, 
we have previously shown that force generation is not a deterministic mechanism but follows a probabilistic 
process and that underlying dynamical events occur on different time scales varying from 100 ms to 5 s 17 . 

For this study we have used optical tweezers to identify the elementary events underlying force generation. 
When an optically trapped bead seals on the lamellipodium membrane, Brownian fluctuations are drastically 
reduced revealing the fine structure of force generation: when a lamellipodium pushes a trapped bead, the 
autocorrelation function p(t) of the bead position decays with multiple time constants up to 50 ms, while during 
Brownian fluctuations p(t) decays with a single time constant less than 1 ms. The distribution of bead velocities 
has long tails with frequent large positive and negative values associated to forward and backward jumps 
occurring in 0.1-0.2 ms. These jumps have varying amplitudes up to 20 nm and their frequency and amplitude 
are reduced when actin turnover is slowed down by the addition of Jasplakinolide 18 and when the action of myosin 

II is inhibited by the addition of Blebbistatin 19 ' 20 . These jumps constitute the elementary events underlying force 
generation. 



Results 

Neurons from dorsal root ganglia (DRG) of P10-P12 rats were isolated and plated on poly-L-lysine-coated glass 
coverslips, positioned on the stage of an inverted microscope used for imaging and force measurement 17 (see 
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Methods). After 24 to 48 hours, lamellipodia emerged from DRG 
soma. Silica beads with a diameter of 1 urn were trapped with an 
infrared (IR) optical tweezer in front of the lamellipodia (Fig. la and 
f): when the lamellipodia protruded and displaced the bead, the 
exerted force F = (F x ,F y ,F z ) was measured with sub pN sensitivity 
at 10 kHz resolution. The bead position x = (x,y,z) was measured 
with a quadrant position detector (QPD) using back focal plane 
(BFP) interferometry 16,21 . Lamellipodia grew by 1 urn within 20- 
30 s and displaced the beads trapped with a low (k x and k y equal to 
0.0155 pN/nm, and k z equal to 0.005 pN/nm) and a high stiffness 
(k x , and k y equal to 0.1 pN/nm and, k z equal to 0.03 pN/nm; Fig. la- 
e). The QPD detects reliably lateral displacements less than 250 nm 
(see Methods) and bead displacements within this range were 
observed with the high trap stiffness. Often lamellipodia pushed 



the bead both laterally and axially (Fig. lf-h) and recordings of 
the bead position became noisier (Fig. Ik). In contrast, when adhe- 
sion forces caused the bead to seal onto the cellular membrane 
of retracting lamellipodia (Fig. li-j) Brownian fluctuations de- 
creased (Fig. Ik). If growth cones were fixed with paraformaldehyde, 
suppressing all cellular motility, no noise increase was observed (see 
Supplementary Methods and Supplementary Fig. SI online). 

During adhesion, the variance (a x 2 , a y 2 ,a z 2 ) could decrease by 5- 
10 times reaching values below 10 nm 2 (Fig. lm) so that the fine 
structure of force generation could be observed. The amplitude 
of the adhesion force F ad was measured as the maximal force before 
the bead returned into the trap 22 " 24 (Fig. 11). Large values of F ad 
reduced more Brownian fluctuations (Fig. In). If F ad is larger than 
30 pN, i.e. the maximal restoring force of the optical trap, when the 
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Figure 1 | During a push, recordings of the bead position become noisier, but not during a pull, (a)-(e) The protruding leading edge of a lamellip odium 
pushes an optically trapped bead by 1 urn within 25 s. (f) A bead trapped in front of a lamellipodium emerging from the soma of a DRG neuron, (g)-(h) 
High resolution images during a push. At 24 s the bead is in the optical trap (g) and when the lamellipodium grows, it pushes the bead (47 s) displacing it 
both laterally and axially (h). (i)-(j) As in (g-h) but during a pull. When the lamellipodium retracted, the bead returned inside the trap (56 s). Following 
bead adhesion, the bead was pulled away from the trap (88 s). Crosses indicate the centre of the optical trap, (k) The three components (x,y,z) of the bead 
displacement. Insets highlight the increase of noise during the push (violet arrows), the decrease of noise during the pull (cyan arrows), and green arrows 
refer to Brownian fluctuations. (1) The three components (x,y,z) of bead displacement during adhesion and retraction in another experiment. At 38 s the 
bead returned into the trap and the adhesion force was measured ( 1 1 pN). (m) Change of variance for the three components in (1). (n) Relation between 
fractional variance reduction and modulus of adhesion force in control conditions (red symbols) and in the presence of 25 nM Jasplakinolide (black 
symbols). The red and black lines represent the linear fit in control conditions and Jasplakinolide, respectively. 
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lamellipodium retracts or vigorously protrudes it could move 
the bead out of the trap. The bonding of a single integrin mole- 
cule to the lamellipodium leading edge is larger than 40 pN 25 and 
therefore, when integrin molecules cause adhesion, the bead is 
strongly attached and will follow the lamellipodium motion also 
out of the optical trap. During adhesion, the power spectrum density 
(PSD) of Brownian fluctuations was fitted by the sum of two or 
three Lorentzian distributions (see Supplementary Methods and 
Supplementary Fig. SI online). 

When actin turnover is reduced no noise increase is observed. 

During a push the variance of bead displacement increased by 2-4 
times (Figs lk and 2a) possibly because of modifications of the 
trap stiffness, fluctuations of adhesion forces or properties of force 
generation. The addition of 100 nM Jasplakinolide, known to reduce 
actin turnover 18 almost completely abolished force generation, 
but a lower concentration of 25 nM slowed down growth cone 
motion without blocking force generation. In the presence of 
25 nM Jasplakinolide, lamellipodia were still able to displace beads 
laterally, but no increase of variance was observed (compare Fig. 2a 
and b). We computed the relation between the lateral component of 
the force (F x or F y ) and the associated changes of variance <j\ 2 . In the 
presence of 25 nM Jasplakinolide, of never increased (grey black 
traces in Fig. 2d; n= 9), but often decreased. On the contrary, in 
control conditions, <5\ increased by 2-4 times (red orange traces 
in Fig. 2d; n= 13). During protrusion the maximal average 
velocity <v max > was 50 nm/s (n = 24), whereas in the presence 
of Jasplakinolide it was 35 nm/s (n =15) 17 . The mean value of the 
modulus of F ad m control conditions was 6.2± 3.3 pN (n=7) and in 
the presence of 25 nM Jasplakinolide was 6.1± 3.1 pN (n=9) 
suggesting that adhesion between the bead and the lamellipodium 
is not affected by Jasplakinolide. Bead displacements and exerted 
forces were very similar in control conditions and in the presence 



of 25 nM Jasplakinolide suggesting that the observed variance 
increase is not caused by local changes of trap stiffness but it is a 
genuine property of force generation. 

When myosin II is inhibited force generation occurs with a 
reduced noise increase. The molecular motor myosin II plays an 
important role in force generation by speeding up actin filament 
disassembly 26 and therefore we analysed the effect of Blebbistatin a 
well known inhibitor of myosin II 19,20 . The addition of 20 uM 
Blebbistatin slowed down lamellipodia motion, but did not abolish 
force generation (Fig. 2c). During lateral push, however, cti 2 did not 
increase by more than 100 % and in some cases (2 out of 6; blue 
traces in Fig. 2d) decreased, but not as observed with Jaslpakinolide 
(Fig. 2d). 

The increase of noise is related to the contact area between the 
bead and the lamellipodium leading edge. The area in direct 
contact A c with a silica bead with a diameter of 1 urn could vary 
from less than 0.1 to up to 1.5 um 2 17 . This contact area mediates all 
mechanical interactions between the bead and force generation 
mechanisms inside lamellipodia. Therefore we have analysed the 
relation between A c and the amplitude of generated force F and 
associated changes of variance <J\. In control conditions, during a 
lateral push it is possible to measure reliably changes of A c and when 
A c increases (Fig. 3a) we have often (7 out of 10 experiments) 
observed a concomitant increase of F (Fig. 3b) and of <J\ 2 (Fig. 3c). 
When the bead went out of the optical trap (broken vertical line) 
the A c - obtained from videomicrographs - remained constant or 
increased while measurements of F and of <j\ 2 are not reliable. 
Therefore, the increase of <j\ 2 observed during lateral pushes is 
caused by the combined effect of force generation and of the asso- 
ciated increase of A c . A c is measured through the objective of the 
microscope, viewing axially the lamellipodium, therefore we were 
not able to determined changes of A c during a vertical push. 
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Figure 2 | Change of noise in control conditions, in the presence of Blebbistatin and in the presence of Jasplakinolide. (a) The longitudinal components 
of the bead displacement during a lateral push in control conditions showing a clear noise increase, (b) As in (a) but in the presence of 25 nM 
Jasplakinolide. No noise increase is observed, (c) As in (a) but in the presence of 20 [iM Blebbistatin. In this case G\ slightly increased in 4 out of 6 cases. In 
2 out of 6 cases <5\ decreased but to a lesser extent in comparison to Jasplakinolide. (d) Relation between force and variance for lateral push in control 
conditions (red shades), in the presence of 20 [iM Blebbistatin (blue shades) and 25 nM Jasplakinolide (black shades). 
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Properties of noise during pushes. Following adhesion, in several 
experiments (n=4 in control conditions; n=6 with Jasplakinolide; 
n=5 with Blebbistatin) the variance of the axial component, a a 2 , 
decreased to less than 6 nm 2 (Fig. 4a) and subsequently, when the 
lamellipodium pushed the bead, a a 2 increased and fluctuations with 
novel properties appeared. 

Visual inspection indicated the existence of rapid discontinuities, 
i.e. of jumps. Therefore, we computed the bead velocity v (Fig. 4b), by 
convolution of bead position with the derivative of a Gaussian func- 
tion, (-t/((2*7t) 1/2 *a 3 ) exp (-t 2 /2*a 2 ) with a value of a varying from 
0.2 to 0.4 ms. During Brownian fluctuations, velocities had a 
Gaussian distribution (Fig. 4c) but not during a push (Fig. 4d): 
indeed their distribution had a central lobe fitted by a Gaussian 
distribution, but had also long tails with large positive and negative 
values. These sudden changes of velocity correspond to forward (j + ) 
and backward (j~) jumps. Similar tails, but less pronounced, could be 
detected also when force generation was not preceded by adhesion. 
During Brownian fluctuations the autocorrelation function p zz (t) of 
bead displacement decayed with a single time constant 0 of 0.64 ± 
0.12 ms (n=20), but duringpushes p zz (t) decayed with multiple time 
constants varying from less than 1 ms up to 50 ms (Fig. 4e). During 
force generation, fluctuations in three coordinates (x,y,z) were more 
correlated: during Brownian fluctuations the cross-correlation p zl (t) 
between z and one lateral component (x or y) decayed with a time 
constant 0 of 0.62 ± 0.15 ms (black traces in Fig. 4g), but during a 
push the value of 0 increased to 6.0 ± 1.4 ms (red traces in Fig. 4g). 
The increase of the time constant 0 of p zz (t) and p z i(t) observed 
during a push was attenuated by 25 nM Jasplakinolide (black trace 
in Fig. 4f and red traces in Fig. 4h). 

All these observations indicate that: i- fluctuations observed dur- 
ing pushes do not originate from thermal motion but are caused by 
the randomness of the elementary events underlying force genera- 
tion; ii - molecular mechanisms underlying force generation are 
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Figure 4 | During pushes the autocorrelation function p(t) of bead position decays with multiple time constants and the distribution of bead velocities 
has long tails, (a) The z component of the bead displacement during Brownian fluctuations (b.m), adhesion and push, (b) Velocity of bead displacement 
in (a), (c) Distribution of velocities during Brownian fluctuations shown in (b). A Gaussian function (red line) fits perfectly the experimental distribution, 
(d) As in (c) during the push shown in (b). (e) Autocorrelation function of vertical bead displacement p zz (t), during Brownian fluctuations, adhesion, and 
push after high pass filtering with a cut-off frequency at 1 Hz (see Supplementary Methods and Supplementary Fig. S2 online). p zz (t) decays with a time 
constant 9 equal to 0.64 ms during b.m. but, during pushing, it has multiple time constants up to 50 ms. (f) The effect of 25 nM Jasplakinolide on p zz (t), 
during pushing (black). The longest time constant of the auto-correlation decreases to 20 ms (red trace), (g) Cross-correlation p zx (t) during b.m. (black 
shades) and during a push (red shades). p zx (t) decays with a time constant 0 equal to 0.62 ms during b.m. and increases to 6.0 ms during a push, (h) The 
effect of 25 nM Jasplakinolide on p zx (t) during pushing. p zx (t) decays with a time constant 0 equal to 0.72 ms during b.m. and increases to 2.47 ms during 
a push. 
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Figure 3 | Concomitant change of variance and contact area during force 
generation, (a) Time evolution of estimated contact area A c (see Ref. 17) 
between the bead and the lamellipodium leading edge, A c , during a push. 
A c at frame i, A c (i), is calculated as A c (i) = 2n [1 — cos(oCi/2)] r 2 , where a { is 
the angle corresponding to the arc of the bead in close contact with the 
leading edge of the lamellipodium and r is the bead radius, as shown in the 
inset, representing a lamellipodium pushing the trapped bead, (b) 
Concomitant time evolution of the force exerted by the lamellipodium 
during the push analyzed in (a), (c) Concomitant time evolution of the 
variance during the push analyzed in (a). The broken vertical line indicates 
the time when the bead is pushed out of the trap. 
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spatially coherent and structured; iii- force generation is character- 
ized by jumps. 

Jumps underlie force generation. Among the vast repertoire of 
algorithms used to detect jumps, we adopted a nonlinear diffusion 
filtering 27 ' 28 , approximating the original data (grey traces in Fig. 5a, b 
and c) with a smooth piece-wise function (red lines) interrupted by j + 
or j" jumps (black vertical segments). This algorithm depends on two 
parameters (see Methods): the contrast X, related to the smallest 
detectable jump, and the scale x determining the temporal window. 
In order to establish our sensitivity and to determine the values of 
X and x we attached a silica bead to the bottom of a coverslip 
which was moved by a piezo manipulator. When the variance of 
displacement fluctuations of the stuck bead was 3.8 nm 2 , as during 
adhesion (Fig. 4a), with the values of 0.5 nm and 0.1 ms for X and x 
respectively, we could detect jumps of 2 nm. With these values of X 
and x, the algorithm detected jumps primarily during pushes (Fig. 5a, 
b and c). 

Immediately after adhesion during the push, (Figs 4a and 5a), we 
detected forward j + and backward j~ jumps ranging from 2 to 20 nm 



(Fig. 5d). Jumps were observed only when lamellipodia pushed the 
bead, but very rarely when beads sealed on the lamellipodia mem- 
brane and retracted, suggesting that jumps do not reflect unspecific 
attachment/detachment events between the lamellipodium and the 
substratum and/or between the actin network and the membrane. 
Jumps lower than 2 nm could not be detected because of noise lim- 
itations. These jumps appear to be the elementary events underlying 
force generation in neuronal lamellipodia. 

Distributions of jumps amplitude in control conditions (Fig. 5d) 
were fitted by the exponential distributions A+ e an( j 
A_ e " C-'j-*) with values of 5.2±1.3 and 4.9±1.2 nm for the mean 
size of positive j + * and negative jumps j~* (n=4). In the presence of 
25 nM Jasplakinolide (Fig. 5c) smaller jumps (Fig. 5f) ranging from 2 
to 8 nm were detected (n=6) and jump distributions were fitted 
by the same exponential distributions but with lower values of j + 
and j~ equal to 2.4 ±0.3 and 2.2 ±0.4 nm, similar in size to the poly- 
merization step size (2.7 nm). If Jasplakinolide is reducing jump 
frequency by stiffening the connection between the lamellipo- 
dium and the bead, it is expected also to modify the adhesion force 
between the bead and the lamellipodium, but this was not observed 
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Figure 5 | Forward and backward jumps are the elementary events underlying force generation, (a) Magnification of the z component of Fig. 4a during 
adhesion and push. Original traces in grey were filtered by the non linear diffusion algorithm (see Methods) providing a smooth component (red curves) 
and jumps (in black). Very few jumps were detected during adhesion but they could be observed very often during a push, (b) Magnification of the z 
component during adhesion and push in the presence of 20 uM Blebbistatin. The original traces (in gray) were filtered as in (a). Jumps with smaller 
amplitude than in control conditions were detected, (c) As in (a) and (b) in the presence of 25 nM Jasplakinolide. In this case jumps with an amplitude 
smaller than the amplitude obtained both in control conditions and in the presence of Blebbistatin were detected, (d)-(f) Density of upward j + and 
downward j~ jumps during push in control conditions (d), in the presence of Blebbistatin (e), and Jasplakinolide (f). These distributions of jump 
amplitude were fitted (black lines in (d) and (e), red lines in (f)) - for values of j + and j~ larger than 2 nm-by the exponential distributions A+e" (i+/j+ * ) and 
A. e " The fitting was performed with the values of of 129 and 128 events/s for the jump frequency of positive and negative jumps, A+ and A_, 

respectively, and 5 and 4.8 nm for the mean size of positive and negative jumps, j + * and respectively (d). In the presence of Blebbistatin the values of 
A+, A_, j + * and j~* were 87 and 80 events/s and 3.5 and 3.3 nm, respectively (e). In the presence of Jasplakinolide the values of A+, A_, j + * and j~* were 
44 and 50 events/s and 2.5 and 2.3 nm, respectively (f). 
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experimentally (Fig. In). Therefore, the effect of Jasplakinolide 
on jump frequency is likely to be caused by a reduced actin turnover. 
In the presence of 20 uM Blebbstatin (Fig. 5b) detected jumps 
have an amplitude ranging from 2 to 15 nm (n=5). Distributions 
of jumps amplitude in the presence of Blebbistatin were fitted by 
exponential distributions with values of j + * and j~* equal to 
3.4 ±0.9 and 3.2 ±0.8 nm. Jumps in the presence of Blebbistatin 
occurred with a frequency 50% lower than in control conditions. 
Therefore Jasplakinolide decreased the jumps frequency and their 
amplitude more than Blebbistatin. 

The detection and measurements of jumps in the presence of noise 
is a difficult (ill-posed) problem 29 that should not be underestimated. 
Therefore, in order to show that jumps are real and not artifacts of 
used algorithms, it is needed to verify that large values of bead velo- 
city v (Fig. 6a) were coincident with jumps (Fig. 6b) detected by the 
non linear diffusion algorithm (see vertical lines in panels 6a and 6b). 
Given the time series of bead position (x n n= 1,..N) the computation 
of the instantaneous velocity does not require any parameter, 
because the velocity is equal to (x n+1 - x n )/At, but two parameters 
are involved in the nonlinear diffusion algorithm (k and x). In order 
to establish co-localization in a quantitative way, large values of v 
were assumed to be those belonging to the tails of the velocity dis- 
tribution outside the Gaussian function fitting its central lobe (see 
Fig. 4d) and these values of v co -localized in a time window At of less 
than 0.3 ms with detected jumps. The analysis of the Rate of True 
Positive (RTP) and of False Positive co-localization 30 (Fig. 6c) indi- 
cates that for At equal to 100 |is, i.e. the used sampling interval, RTP 
is larger than 80% and becomes close to 100% for At equal to 300 jas. 
This analysis indicates that jumps detected by the non linear dif- 
fusion algorithm co-localize very precisely with large values of bead 
velocity. 

Jumps were clearly detected when force generation developed fol- 
lowing adhesion, i.e. when Brownian fluctuations were reduced 
(Fig. 5). However, more often force generation developed without 
being preceded by bead adhesion (Fig. 7a). Therefore we asked 
whether it was possible to determine the existence of jumps also 
when force generation did not follow bead adhesion. As force 



generation is characterized by a large value of the autocorrelation 
function (Fig. 4) we computed for all three components x,y and z, 
Pxx(t)> Pyy(t), p Z z(t). During Brownian fluctuations p^t) are expo- 
nentially distributed with a value of x less than 1 ms, but during push 
pii(t) become broader decaying with several time constants (Fig. 4). 
Therefore we computed the integral Ci(t) of pii(t) at each time and 
force generation was identified to occur when Q(t) increased by at 
least 10 times (Fig. 7b). Under these circumstances the variance of 
bead fluctuations at the peak of force generation was significantly 
higher than during Brownian fluctuations (Fig. 7c). During this 
phase the central lobe of the distributions of dx/dt, dy/dt and dz/dt 
was fitted by a Gaussian function, but tails corresponding to large 
positive and negative velocities were detected (see arrows in Fig. 7d- 
f) indicating the existence of forward and backward jumps. Jumps 
detected by the non linear diffusion algorithm during these events 
(Fig. 7g-i) have amplitudes ranging up to 20 nm, as those detected 
after adhesion (Fig. 5). 

If jumps are the elementary events underlying force generation 
their sum must be close to the observed net protrusion. Therefore we 
compared the net protrusion Prot(At) in the time window At with 
the sum of all forward jumps (L At j + ) minus the sum of all backward 
jumps (Z At j ) occurring in At: P(At) (red line) was very similar to 
2 At j + - 2 At j~ (black line in Fig. 8a and b) in control conditions as 
well as in the presence of Jasplakinolide. We next asked whether force 
generation developed by an increase of the frequency of jumps, i.e. 
A + and A_ or by their mean amplitude, i.e. j + * and j"*. Therefore we 
estimated A+, A_, j + * and j~* in 0.5 s intervals during force genera- 
tion: force generation developed by a combination of an increase of 
jumps frequency and of their mean values. The observation that 
forward and backward jumps sum to net protrusion is an additional 
test for internal consistency of the used procedure for jump detec- 
tion, providing further support that jumps are the elementary events 
underlying force generation. 

Discussion 

The results of the present manuscript show that force genera- 
tion in neuronal lamellipodia of rat DRG neurons, is composed by 
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Figure 6 | Colocalization of jumps and large values of bead velocity, (a) Bead velocity during a push obtained by the convolution of the bead 
displacement with the derivative of a Gaussian function (— t/(2*7i) 1/2 * a 2 ) exp (— t 2 /2*a 2 ) with a = 0.1 ms. A jump and a large value ofvcolocalizeifthey 
occur in a time window At of less than 0.3 ms. Large values of v were those belonging to the long tails of velocity distribution. These velocities had an 
absolute value larger than 3 times the standard deviation of the Gaussian fitting the central lobe of the velocity distribution (see Fig. 4d). (b) Jumps 
detected by nonlinear diffusion in the same portion of the push shown in (a), where velocity was computed. Original trace, gray; smoothed component, 
red; detected jumps, black. Red and black dotted lines highlight colocalization of positive j + and negative j~ jumps, respectively with large values of v. 
(c) Rate of True Positive colocalization and of False Positive colocalization for increasing values of At from 0 to 1 ms, Asterisks represent the classifiers in 
which positive (red) and negative (black) jumps colocalize exactly with large values of v 30 . 
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Figure 7 | Characterization of force generation during a push in the absence of adhesion. Episodes of force generation were identified as increases of at 
least 10 times of the integral Q(t) of pii(t), i=x,y,z. (a) The three components (x,y,z) of the bead displacement during a push, (b) Integral C x (t) (blue), 
C y (t) (green), and C z (t)(red) of the autocorrelation function pxx(t), pyy(t), and p zz (t) of each of the three components of the bead displacements shown in 
(a), (c) Change in time of the bead displacement variance for the three components in (a). Variance computed in time windows of 0.1 s after high pass 
filtering at 1 Hz. (d)-(f) Distribution of values of velocity dz/dt (d), dx/dt (e), and dy/dt (f), during force generation shown in (a). The black line 
represents the Gaussian fit to the distribution. The arrows highlight the tails associated to forward and backward jumps, (g)-(i) Density of upward j + (red 
histograms) and downward j~ (blue histograms) jumps during the push shown in (a) for the z (g), x (h), and y (i) component, respectively. 




5 s 5 s 



Figure 8 | The sum of forward and backward jumps is equal to the net protrusion, (a) Bead vertical displacement during a push (red line) in control 
conditions. The black line represents the sum of all forward jumps (Z At j + ) minus the sum of all backward jumps (Z At j") occurring in the time window 
At =0.5 s, calculated over the whole push, (b) As in (a) in the presence of 25 nM Jasplakinolide. 
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elementary events corresponding to forward and backward jumps. 
These jumps have an amplitude ranging from 2 to 20 nm and have 
varying orientation in the 3D space. These jumps are not observed 
when growth cones were fixed with paraformaldehyde, suppressing 
all cellular motility (Supplementary Methods and Supplementary 
Fig. SI online) and their amplitude and frequency were reduced by 
treating growth cones with 25 nM Jasplakinolide (Fig. 5). Jumps 
detected by the nonlinear diffusion algorithm colocalize with high 
values of the instantaneous bead velocity (Fig. 6) and the net protru- 
sion of lamellipodia is the net sum of forward and backward jumps 
(Fig. 8). For all these reasons, jumps - here described-are neither 
artifacts of the detection procedure nor are caused by changes of 
properties of the optical trap. Detected jumps represent the element- 
ary events underlying force generation in DRG lamellipodia. 

Force generation occurs at different rates. At the slowest rate the 
lamellipodium leading edge advances smoothly with forward and 
backward jumps (Fig. 5f) with an amplitude similar in size to the 
mean polymerization step size (2.7 nm) observed during actin fila- 
ment polymerization 31 ' 32 suggesting that actin monomers are added 
one by one to the existing network of filaments. At the fastest rate 
(Fig. 5d) larger jumps are observed and they could result from the 
insertion of small actin oligomers 33 and by the occurrence of a burst 
of actin polymerization in a single or neighboring actin filaments. 
Oligomers of 5-10 actin molecules are present in lamellipodia, prim- 
arily as a result of actin filaments depolymerization and could be used 
for actin filament assembly, as in yeast Saccharomyes cervisiae 33 . 
Experimental determinations of free G-actin in lamellipodia varied 
from 1-3 uM 9 " 12 to values 100 higher 34 . In vitro determination of 
actin polymerization rates provides values ranging from 11.6 to 
38 (1/uM s) 912 , but the bulk turnover of actin subunits is 100-200 
times faster in cells than in vitro 1,35,36 . Therefore, the occurrence of a 
fast and vigorous polymerization rate of a single actin filament is 
possible. As the estimated density of actin filaments impinging upon 
the leading edge of a lamellipodium is between 100 and 200 per 
urn 2 13,32 a burst of polymerization of several actin filaments and 
an appropriate spatial environment could also produce a discrete 
forward step of 10-20 nm of the lamellipodium leading edge. 
When actin turnover is reduced by Jasplakinolide 18 and when myosin 
II is inhibited by Blebbistatin 19,20 , force generation still occurs but at a 
slower rate (Fig. 5e and f). Jumps produced by DRG lamellipodia 
do not have a discrete amplitude as observed in the actin-based 
movement of the bacterium Listeria monocytogenes, with a predom- 
inant jump size of 5.4 nm 37 . 

Force generation in lamellipodia depends on several factors such 
as the availability of actin monomers/oligomers, the presence of 
molecular motors such as myosin II and a large variety of controlling 
proteins 6 . This complexity is at the basis of the observed dynamics, 
reminiscent of self organized criticality 38 . 

Methods 

Neuron preparation. Wistar rats (P10-P12) were anaesthetized with C02 and 
sacrificed by decapitation in accordance with the Italian Animal Welfare Act. The 
Ethics Committee of the International School for Advanced Studies (SISSA-ISAS) has 
approved the protocol (Prot.n. 2189-II/7). Dorsal Root Ganglia (DRG) were 
incubated with trypsin (0.5 mg/ml, Sigma- Aldrich, Milan, Italy), collagenase (1 mg/ 
ml, Sigma- Aldrich), and DNase (0.1 mg/ml, Sigma-Aldrich) in 5 ml Neurobasal 
medium (Gibco, Invitrogen, Milan, Italy) in a shaking bath (37°C, 35-40 minutes). 
DRGs were mechanically dissociated, centrifuged at 300 rpm, resuspended in culture 
medium and plated on poly-L-lysine-coated (0.5 /ig/ml, Sigma-Aldrich) coverslips. 
Cells were incubated for 24 to 48 hours followed by the addition of nerve growth 
factor (50 ng/ml; Alomone, Israel) before the measurements. 

Optical tweezer set-up. The optical tweezers set-up was built as previously 
described 17 ' 39 . The dish containing the differentiating neurons and the beads (PSI- 
1.0NH2, G.Kisker GbR, Steinfurt Germany) was placed on the microscope stage 
which could be moved by a 3 axis piezoelectric nanocube (17 MAX 301, Melles Griot 
Inc., USA). The temperature of the dish was kept at 37°C by a Peltier device. The bead 
position was determined in the x,y and z plane with a lateral and axial accuracy of 2 
and 5 nm respectively, which was obtained from the analysis of the interference 
between forward scattered light from the bead and unscattered light 16 ' 21 . The back 



focal plane of the condenser was imaged onto a QPD (C5460SPL 6041, Hamamatsu, 
Milan, Italy) and the light was converted to differential outputs digitized at 10 kHz 
and low pass filtered at 5 kHz. Both the lateral and axial trap stiffness, kxy = (kx, ky) 
and kz, respectively, as well as the detector sensitivity were calibrated using the power 
spectrum method 16 with voltage signals filtered and digitized at 5 kHz. In order to 
reduce and possibly avoid all mechanical perturbations affecting the measurement of 
x= (x,y, z), the optical tweezers set-up was kept in an isolated and sound-proof room 
and the scientists performing the experiments, controlled all operations remotely 
from a separate room. In order to have good mechanical stability it was necessary to 
position all power supplies of used equipment in a separate room and to avoid flying 
cables by properly securing them. In this way we reduced perturbations, which could 
have affected previous investigations. 

Jumps determination by non linear diffusion filtering. In order to detect jumps, we 
used an algorithm based on non linear diffusion 27 ' 28 . After selecting the part of the 
trace of interest, the original signal was approximated with a smooth piece-wise 
function where the discontinuities, i.e. rapid and large changes, were identified as 
jumps. The non-linear diffusion is an iterative process based on the choice of two 
parameters: the contrast X, related to the minimal jump amplitude detected, and the 
scale x, determining the temporal window of jumps. The values of X and x were set 
equal to 0.5 nm and 0.1 ms, respectively, so to detect 2 ground truth jumps (see 
Supplementary Methods and Supplementary Fig. S3 online). The algorithm is based 
on the Toolbox of Frederico D'Almeida (see http://www.mathworks.com/ 
matlabcentral/fileexchange/3710-nonlinear-diffusiontoolbox). 

We compared the detection of jumps using the same values of X and x from traces 
obtained in different conditions (Brownian fluctuation recording, adhesion, and 
push; see Supplementary Methods and Supplementary Fig. S3 online), but scaled so to 
have the same width of the central lobe of the velocity distribution. During pushes, 
jumps were detected with a rate about 4 times higher than in the other conditions. 
Therefore, if the variance of Brownian fluctuations of the trapped bead decreases to or 
below 4 nm 2 , the overall system can detect reliably 2 nm jumps. 
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